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We present a direct measurement of the density profile of a two-dimensional Mott Insulator 
formed by ultracold atoms in an optical lattice. High resolution absorption imaging is used to 
probe the "wedding-cake" structure of a trapped gas as it crosses the boundary from a unit-filled 
Mott insulating phase to the superfluid phase at finite temperature. Detailed analysis of images 
yields measurements of temperature and local compressibility; for the latter we observe a strong 
suppression deep in the Mott-insulating phase, which is recovered for the superfluid and normal 
phases. Furthermore, we measure spatially resolved fluctuations in the local density, showing a 
suppression of fluctuations in the insulator. Results are consistent with the fluctuation-dissipation 
theorem for insulator, superfluid and normal gas. 
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The observation of the Superfluid (SF) to Mott- 
insulator (MI) phase transition of ultracold atoms in op- 
tical lattices [1] marked an experimental breakthrough 
in the study of many-body physics in an atomic sys- 
tem. This development provided the first tangible ex- 
ample of a quantum phase transition (one that occurs 
even at zero temperature) and suggested a highly cor- 
related and gapped Mott insulating state of a Bose gas, 
which persists down to zero temperature without Bose- 
condensation and macroscopic phase coherence. Since its 
theoretical inception [H, H, 01 , two of the most celebrated 
properties of the bosonic Mott insulator have been its 
incompressibility and suppression of local density fluc- 
tuations [5], induced by enhanced inter-particle interac- 
tion and reduced mobility in optical lattices. The result 
for a trapped atom gas, where the local chemical po- 
tential varies in space, is the remarkable "wedding-cake" 
density profile, where successive MI domains are mani- 
fest as plateaus of constant density. Relatedphenomena 
have been studied through the coherence trans por t 
P, Q, noise correlations [8], and number variance [9l,]lQt, 
but direct observation of the incompressibility has proven 
difficult due to the inhomogeneous nature of all experi- 
ments to date, and to the technical difficulty of spatially 
resolved measurements. Innovative experimental efforts 
have yielded evidence [ll|, [l2|] that these plateaus exist, 
though none has directly observed this effect by imaging 
a single physical system in situ. 

We report studies based on direct in-situ imaging of an 
atomic MI. By loading a degenerate Bose gas of cesium- 
133 atoms into a thin layer of a two-dimensional optical 
lattice potential, and adiabatically increasing the optical 
lattice depth, we observe the emergence of an extremely 
flat density near the center of the cloud, which corre- 
sponds to a MI phase with accurately one atom-per-site. 
From density profiles, we extract important thermody- 
namic and statistical information, confirming the incom- 
pressibility and reduction of density fluctuations in the 
MI as described by the fluctuation-dissipation theorem. 



The single layer, two-dimensional (2D) optical lattice is 
formed by two pairs of counter-propagating laser beams 
derived from a Yb fiber laser at wavelength A = 1064 nm. 
The pairs are oriented orthogonally on the horizontal 
{x — y) plane, forming a square lattice with site spac- 
ing d = \/2 = 0.532 jiia. A weak harmonic potential of 
Vh = m{uo1x^ -^uOyy'^)/2 localizes the sample, where m 
is the cesium mass, and the geometric mean of the trap 
frequencies is uOr = ^J^x^y = 27r x 9.5(1 + y/82^i^)Hz; 
we have included a weak dependence on the lattice depth 
V {Er = h 1.3kHz is the lattice recoil energy and h is 
Planck's constant.) Vertical confinement is provided by 
an additional vertical optical lattice with a site spacing 
4 /im, formed by two beams intersecting at an angle of 
15°, confining atoms to an oscillator length = 0.30 jam. 
The sample is loaded into a single site of the vertical lat- 
tice, kept deep to prevent vertical tunneling. Tunneling 
in the horizontal 2D lattice is controlled by varying the 
lattice depth V f]\ . Details on preparation of the atomic 
sample can be found in the supporting material and Ref. 

0. 

We obtain a top view of the sample using absorption 
imaging, directly revealing the atomic surface density 
n(x, y) on the horizontal plane. The imaging resolution is 
3 ~ 4/im, and magnification such that one imaging pixel 
corresponds to an area of (2/im)^ on the object plane. 
Conveniently, unit filling in a 2D optical lattice has a 
moderate optical density of OD = 0(1) on resonance. 

The superflulid-to-Mott insulator (SF-MI) transition 
of ultracold atoms in an optical lattice is described by 
the Bose-Hubbard model, characterized by on-site inter- 
action U and the tunneling J 0. In 2D optical lattices, 
superfluid is converted into a MI when U/J exceeds 16 
@, Here, the SF-MI phase transition can be induced 
by either increasing the lattice potential V P, 0, [l2[ or 
the atomic scattering length a via a magnetically-tuned 
Feshbach resonance [15], together providing complete, in- 
dependent control of U and J. 

Atomic density profiles in the lattice are shown in 
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FIG. 1: False color density profiles and line cuts oi N = 7500 
ultracold cesium atoms at scattering length a = 310 as in 2D 
optical lattices. (A) Superfluid regime (shallow lattice V = 
0.14^k), (B) Phase transition regime (medium lattice depth 

V = 7.2Er), and (C) Mott insulator regime (deep lattice 

V = 3SEr) The red dashed line indicates the expected density 
for one atom-per-site. Line cuts pass through the center of 
the sample. 



Fig. [TJ For weak lattice depths (superfluid regime), the 
density profiles are bell-shape, with negative curvature 
at the center (Fig. [T]A.), indicating a finite, positive com- 
pressibility dictated by the interaction coupling constant 
(discussed below.) In sufficiently deep lattices, we ob- 
serve a flattened density at the center of the sample 
(Fig. [Tp, C), indicating development of a Mott insu- 
lating phase with one particle per lattice site. This den- 
sity plateau, an important feature of the MI phase, arises 
from incompressibility. 

A primary check on the MI is to compare the measured 
density in the plateau to that corresponding to one atom- 
per-site, given by MI physics as a "standard candle" of 
atomic density. Using the known scattering cross-section, 
correcting for saturation effects (see Methods), we de- 
termine the plateau density to be n = 3.5(3)//im^, in 
agreement with the expected value = 3.53//im^. 

To distinguish a MI from superfluid or normal gas, we 
histogram the occurrence of pixels h{n) in the images 
corresponding to a density n with a bin size of An <C n. 
The MI plateau, containing a large number of pixels with 
similar atomic density, appears as a peak at n = 
(Fig.[2]A). In general, the occurrence of a particular den- 
sity n can be regarded as the rate at which local chemical 
potential changes with density, multiplied by the number 
of pixels w{/jj)A/j. corresponding to a chemical potential 
between fi and /i + A/i. The occurrence at density n 
is then h{n) = Anw{jii)Ajii/ An ^ Anw{ii)n~^ , where 
n = dn/dji is the local compressibility In a har- 

monic trap, w{ii) = 2t: /mdPuo'^ is constant, and the his- 
togram is a particularly useful tool to distinguish differ- 
ent phases. For a pure BEG in the Thomas- Fermi limit, 
the compressibility is constant to the maximally- allowed 
density n^/c, and results in a constant h{n) for n < ripk 
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FIG. 2: Histograms of density profiles in the MI regime (A, 
V = 38Er, a = 4:60a b) and the superfluid regime (B, 1/ = 
0.5Er, a = 460aB.) The histograms are based on an average 
of three density images. The bin size is An = 0.03. 



(see Fig. EB for < n< For the MI, the 

density is insensitive to chemical potential in a narrow 
range near n = indicating a vanishing compress- 

ibility, and thus a sharp histogram peak at n = l/cP. 
The peak's presence in Fig. [2]A is thus directly related 
to the incompressibility in the Mott phase. Finally, the 
compressibility of a normal (ideal) gas is proportional to 
its density, thus h{n) ex 1/n, leading to the strong upturn 
at low densities in Fig. [2]A,B for both regimes. 

Much more information can be obtained from the den- 
sity profiles, as recently suggested in Ref. [17]. For ex- 
ample, the compressibility in a two-dimensional cylin- 
drically symmetric trap can be written hi = dn/dji = 
—n^{r)/{rmu'^), where we have assumed the local density 
approximation, and that the chemical potential depends 
on the trapping potential ji = jiQ — Vf/(r). For a BEG in 
the Thomas- Fermi regime, the compressibility is a posi- 
tive and constant kbec = 1/^, where g = \^ah'^ /maz 
is the (2D) interaction parameter [lH|. We can thereby 
relate the measured compressibility to that of a BEG as 



f^BEC 



rd ^ az hujr 



(1) 



We evaluate from azimuthally averaged density pro- 
files (Fig. [3] A). Eccentricity of the trap is corrected by 
rescaling the principal axes. Due to the singular nature 
of n'{r)/r near the center, we evaluate n there by fit- 
ting n(r) to a quadratic, n(r) = n(0) — ar^ . The cur- 
vature then gives the compressibility as niO) = 2a/ma;^, 
for which we obtain n/KsEC = 0-4(1) in a weak lattice 
and n/nBEC = 0.006(6) in a strong lattice (See Fig. [3]). 
In the weak lattice (SF regime), the finite and constant 
compressibility at the center agrees with expectation for 
the superfluid phase, though lower than expected, which 
we attribute to finite temperature and calibration of trap 
parameters. The finite temperature is also clear in the 
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FIG. 3: (A) Radially averaged profiles (3 images) in the su- 
perfluid (blue open circles: V — 0.25Er) and MI (red circles: 
V = SSEr) regimes. Particle number and scattering length 
are N = 8100 and a = 460aB, respectively. A quadratic fit to 
the center portion of the sample extracts the curvature near 
r = 0. (B) Normalized compressibilities derived from (A) us- 
ing Eq. (1) in the superfluid (blue open circles) and MI (red 
circles) regimes. The horizontal lines indicate the compress- 
ibility near r = 0, estimated from the quadratic fits in (A). 
The rising compressibility at r = 40d marks the MI boundary. 
The inset shows the dependence of compressibility on atomic 
density. The linear dependence at low densities is fit by the 
solid lines. 



FIG. 4: The fluctuation of local density extracted from a set 
of twelve absorption images in the deep (MI) and weak (SF) 
lattice regimes. The insulator and superfluid show a pro- 
nounced difference at the density of one atom-per-site, where 
the insulator's fluctuation is suppressed by its incompressibil- 
ity. In the superfluid phase, constant compressibility initiates 
a flattening of fluctuations with density. At low densities, in 
both regimes, the fluctuation shows a characteristic ^/n de- 
pendence. The dashed line shows best fit a/ti dependence for 
the normal gas. The total number of atoms was N = 8300 
(SF) and N = 9600 (MI) with a = SlOas for both sets. 



exponential tail of the density profile and the compress- 
ibility, from which we derive the temperature 11(2) nK 
in the superfluid regime and 24(4) nK in the MI regime. 

In a deep lattice (MI regime), we observe a strong re- 
duction of the compressibility in the trap center, below 
that in the superfluid phase for the weak lattice, strongly 
supporting the emergence of a MI phase at the center of 
the sample. Away from center, suddenly increases at 
r = 25(i, then decreases for r > 50d. The exponential 
decay is again consistent with a normal gas. Between MI 
and normal gas {25d < r < bOd), a more detailed mea- 
surement and model of compressibility would be neces- 
sary to identify the local phase. 

Incompressibility necessarily implies a low particle 
number fluctuation; this relationship, a result of the 
fluctuation-dissipation theorem (FDT) (see e.g. [l9, 20]), 
takes the form 

5n^ = KkBT (2) 

Resolved in-situ imaging provides an enticing opportu- 
nity to measure fluctuations of the local density [5, 21], 
and thus check the validity of the FDT. We measured 
fluctuations by recording multiple absorption images, cal- 
culating the variance of density measured in each pixel 
(each collects signals from a patch of (2/im/(i)^ ^ 14 lat- 
tice sites). Fig. m shows the recorded fluctuations, where 



pixels are binned according to their mean density. Fluc- 
tuations consist of detection (photo-electron shot) noise 
and thermal and quantum atomic density fluctuations. 
Detection shot noise can be well-calibrated and modeled 
by analyzing portions of the images with low density; ex- 
tension to higher optical depth (density) shows the weak 
dependence illustrated (Fig. 2]). 

Above the detection noise, density fluctuations (see 
Fig. [4]) show a strong qualitative agreement with the 
compressibility presented in Figure [3] as expected from 
the FDT. For example, the Mott-phase shows a strong 
suppression of fluctuations at the density of one atom- 
per-site. The superfluid regime lacks this feature, instead 
showing a pronounced flattening as the sample transi- 
tions from normal gas to superfluid, as expected from the 
constant compressibility in the superfluid phase (Figure 
[3l inset). Finally, at low density, the normal gas shows 
a temperature-independent fluctuation of 5n = ^\/n^ 
which can be anticipated from the inset of Figure O and 
agrees with the FDT. The coefficient 7 is roughly con- 
sistent with the FDT, and measured imaging resolution 
(see methods). 

Clearly, in situ imaging of the Mott insulator is a pow- 
erful new tool to investigate new quantum phases of cold 
atoms in optical lattices. From the density profiles, not 
only can one observe the density plateau, incompressibil- 
ity and reduction of fluctuations in the Mott insulating 
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phase, but also demonstrate a qualitative validation of 
the fluctuation-dissipation theorem. Relatively modest 
extension of this work holds new promise for studying 
the role of quantum fluctuations, correlation and ther- 
modynamics near a quantum phase transition. 

We thank T.L. Ho, R. Scalettar, E. Mueller, and R. 



Hulet for helpful discussion. This work is supported by 
NSF No. PHY-0747907, NSF-MRSEC DMR-0213745, 
and ARO No. W911NF0710576 with funds from the 
DARPA OLE Program. N.G. acknowledges support from 
the Grainger Foundation. 



[1] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and 

I. Bloch, Nature 415, 39 (2002). 
[2] M. I. Kaganov and A. V. Chubukov, Uspekhi Fizich- 

eskikh Nauk 153, 537 (1987). 
[3] M. P. A. Fisher, P. B. Weichman, G. Grinstein, and 

D. Fisher, Physical Review B 40, 546 (1989). 
[4] D. Jaksch, C. Bruder, J. 1. Cirac, C. W. Gardiner, and 

P. Zoller, Physical Review Letters 81, 3108 (1998). 
[5] B. Capogrosso-Sansone, N. V. Prokof 'ev, and B. V. Svis- 

tunov. Physical Review A 75, 013619 (2007). 
[6] L B. Spielman, W. D. Phillips, and J. V. Porto, Phys 

Rev Lett 100, 120402 (2008). 
[7] M. Kohl, H. Moritz, T. Stoferle, C. Schori, and 

T. Esslinger, Journal of Low Temperature Physics 138, 

635 (2005). 

[8] S. Foiling, F. Gerbier, A. Widera, O. Mandel, T. Gericke, 

and L Bloch, Nature 434, 481 (2005). 
[9] M. Greiner, O. Mandel, T. W. Hansch, and L Bloch, 

Nature 419, 51 (2002). 
[10] F. Gerbier, S. Foiling, A. Widera, O. Mandel, and 

L Bloch, Physical Review Letters 96, 090401 (2006). 
[11] S. Foiling, A. Widera, T. Muller, F. Gerbier, and L Bloch, 

Phys Rev Lett 97, 060403 (2006). 
[12] G. K. Campbell, J. Mun, M. Boyd, P. Medley, A. E. 

Leanhardt, L. G. Marcassa, D. E. Pritchard, and W. Ket- 

terle. Science 313, 649 (2006). 
[13] C. L. Hung, X. B. Zhang, N. Gemelke, and C. Chin, 

Physical Review A 78, 011604 (2008). 
[14] L B. Spielman, W. D. Phillips, and J. V. Porto, Physical 

Review Letters 98, 080404 (2007). 
[15] C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, 

arXiv.org:0812.1496 (2008). 
[16] M. Rigol, G. G. Batrouni, V. Rousseau, and R. T. Scalet- 
tar, arXiv:0811.2219 (2008). 
[17] T.-L. Ho and Q. Zhou, arXiv:0901.0018vl (2008). 
[18] B. Tanatar, A. Minguzzi, P. Vignolo, and M. P. Tosi, 

Physics Letters A 302, PH S0375 (2002). 
[19] G. G. Batrouni, V. Rousseau, R. T. Scalettar, M. Rigol, 

A. Muramatsu, P. J. H. Denteneer, and M. Troyer, Phys- 
ical Review Letters 89, 117203 (2002). 
[20] K. Huang, Statistical mechanics (Wiley, New York, 

1963). 

[21] J. Esteve, J. B. Trebbia, T. Schumm, A. Aspect, C. L 
West brook, and L Bouchoule, Physical Review Letters 
96, 090401 (2006). 



5 



Supporting Material 

Preparation of BEC in a thin 2D optical lattice 

The ^^^Cs BEC is formed in a crossed-beam dipole trap 
by an efficient evaporative cooling method [1]. The dipole 
trap consists of three beams on the horizontal plane: two 
orthogonal beams at the wavelength of 1064 nm (Yb fiber 
laser, YLR-20-1064-LP-SF, IPG), focused to l/e^ radii 
of 350 /im, and one CO2 laser beam (not shown in Fig 
[5j) at the wavelength of 10.6 fim (Gem- Select 100, Co- 
herent), focused to a vertical 1/e^ radius of 70 /im and 
horizontal of 2 mm (see Fig. [5]). The CO2 beam inter- 
sects the Yb laser beams at an angle of 45° and provides 
an enhanced vertical confinement to support the atoms 
against gravity. With = 10^ atoms in a pure con- 
densate, the Thomas-Fermi radii of the condensate are 
{rx.ry.Tz) = (23,14,3.6)/im. 




FIG. 5: Experimental apparatus used to investigate quantum 
phases in two-dimensional optical lattices. Two pairs of coun- 
terpropagating beams (red) form a 2D square optical lattice, 
into which a cesium Bose-condensate is loaded. Atoms are 
confined in the vertical direction by loading into a single site 
of a separate one- dimensional lattice in the vertical direction, 
formed by two beams (blue) intersecting at a small angle. 
A high numerical aperture imaging lens images the shadow 
(absorption) cast by the atoms illuminated by imaging light 
propagating vertically. 

After a pure BEC is obtained, the sample is com- 
pressed vertically by introducing a vertical lattice, 
formed by two laser beams (Innolight Mephisto) inclined 
at +7.5° and -7.5° relative to the horizontal plane. The 
vertical lattice has a spacing of 4 /im and, together with 
the crossed dipole trap, forms an array of 2D oblate "pan- 
cake" potentials. 

In order to load the condensate into a single pancake 
trap, we first ramp the magnetic field to 17.2 G in 400 ms, 
reducing the s-wave scattering length to a <10 a^, and 
then turn on the vertical lattice in 100 ms. Atomic pop- 
ulation in other lattice sites, if any, can be identified by 



observing an interference pattern in time-of-fiight images 
taken from the side. For this work, we observe a suf- 
ficiently weak interference pattern contrast to conclude 
> 98% of the atoms are in a single pancake trap. After 
the vertical lattice is fully turned on, the CO2 laser in- 
tensity is ramped to zero in 100 ms and the scattering 
length ramped to a desired value by tuning the external 
magnetic field. 

The 2D lattice potential in the horizontal {x- and y-) 
directions is formed by introducing retro-reflections of the 
1064 nm dipole trap beams. A continuous evolution from 
a pure dipole trap (with zero retro-reflection) to a 2D op- 
tical lattice (with significant retro-reflection) is achieved 
by passing each dipole trap beam (after it passes through 
the atomic cloud once) through two acousto-optic mod- 
ulators (AOMs) controlled by the same radio-frequency 
(rf) source, then off a retroreflection mirror. The AOMs 
induce an overall zero frequency shift, but permit a dy- 
namic control of the retroreflection intensity over six or- 
ders of magnitude. Onsite interaction energy U and tun- 
neling rate J are evaluated from the measurements of 
the lattice vibration frequencies and the band structure 
calculation. 



Calibration of atomic surface density 

By varying the intensity of the imaging beam, we 
measure the optical depth in the density plateau using 
OD = ln(Mo/M), where M is the number of photons 
collected by a CCD pixel in the presence of the atoms 
and Mo is that without the atoms. The optical depth 
in the plateau is extracted from a fit to the peak in 
the histogram. We then fit the variation of peak optical 
depth assuming OD = na/{l + Mo/Mgat) to determine 
the depth in the zero intensity limit Mq 0, and thus 
the surface density of the sample. Here, cr=0.347 /im^ 
is the known cesium atom-photon cross-section while the 
fit parameter Mgat represents the photon number on a 
CCD pixel at the atomic saturation intensity. 



Fluctuation of atomic density 

The fluctuations in the absorption images are esti- 
mated by taking the average of 11 images under the same 
experimental procedure, and calculating the mean and 
variance of optical depth measured at each CCD pixel. 
Fluctuations are presumed to arise from optical shot 
noise, thermal atomic fluctuation, and long lengthscale 
variations arising from total atom number fluctuation. 
The optical shot noise is calibrated by examining regions 
with negligible atomic density , and ex tended to higher 
optical depth using SODqs cx Vl + e^^ . For the thermal 
cloud, with density n < 0.3 atoms/site, the fluctuation- 
dissipation theorem predicts SNa = with Na the 
number of atoms measured in a given region. This re- 
sult should be valid for a region significantly larger than 



the correlation length, which we expect for the normal 
gas to be on order of the deBroglie thermal wavelength, 
expected to be < 1.5/im for our sample. Though each 
imaging pixel corresponds to an area in the object plane 
consisting of 14 sites, imperfect imaging resolution is 
expected to effectively average away a certain fraction 
of the total fluctuation. This effect can be calculated, 
assuming statistical independence for each site, by sum- 
ming the weight Wij of a resolution-limited spot falling 
within a given pixel j for each lattice site i, giving a vari- 
ance reduced by wf j . The result for our parameters is 
a reduction to 5n = j^/n^ with 7 ~ 0.11(1). This should 
be compared with the fraction of the total fluctuation 
shown in Fig. [H corresponding to thermal fluctuations 
in the superfluid regime. To make this comparison, we 
reject global fluctuations associated with variation of the 
total atom number by subtracting the variance we cal- 
culate after first applying a resolution-spoiling gaussian 
blur to the images from the variance without modifica- 
tion. We find, for the remaining high spatial-frequency 
fluctuations, a best fit to 7 of 0.15(2), using a gaussian 
blur 1/e^ radius of = 14/im to remove global varia- 
tions (the result varies within stated error for blur radii 
7/im< Tb < 28/im). The remaining discrepancy is likely 
due to calibration of imaging resolution, and possibly the 
effect of a nonnegligible correlation length. 
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